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The crossover from weak to strong coupling for a three dimensional continuum model of fermions 
interacting via an attractive contact potential is studied above the superconducting critical temper- 
ature T c . The pair- fluctuation propagator, the one- loop self-energy, and the spectral function are 
investigated in a systematic way from the superconducting fluctuation regime (weak coupling) to 
the bosonic regime (strong coupling). Analytic and numerical results are reported. In the strong- 
coupling regime, where the pair fluctuation propagator has bosonic character, two quite different 
peaks appear in the spectral function at a given wave vector, a broad one at negative frequencies and 
a narrow one at positive frequencies. The broad peak of the spectral function at negative frequen- 
cies is asymmetric about its maximum, with its spectral weight decreasing by increasing coupling 
and temperature. In this regime, two crossover temperatures T* (at which the two peaks in the 
spectral function merge in one peak) and Tq (at which the maximum of the lower peak crosses zero 
frequency) can be identified, with T c <C T * < T*. By decreasing coupling, the two-peak struc- 
ture evolves smoothly. In the weak-coupling regime, where the fluctuation propagator has diffusive 
Ginzburg-Landau character, the overall line-shape of the spectral function is more symmetric and 
the two crossover temperatures approach T c . The systematic analysis of the spectral function iden- 
tifies specific features which allow one to distinguish by ARPES whether a system is in the weak- or 
strong-coupling regime. Connection of the results of our analysis with the phenomenology of cuprate 
superconductors is also attempted and rests on the recently introduced two-gap model, according to 
which a crossover from weak to strong coupling is realized when moving in the Brillouin zone away 
from the nodal points toward the M points where the d-wave gap acquires its maximum value. 



PACS numbers: 74.20.-Z, 74.25.-q, 74.25.Jb 
I. INTRODUCTION 

High-T c cuprate superconductors are characterized by 
doping- and temperature-dependent anomalous proper- 
ties in the metallic and superconducting phases. At low 
doping (i.e., in the underdoped region of the phase di- 
agram), the cuprates display a pseudogap in the single- 
particle excitation spectra and in the spin susceptibil- 
ity, above the superconducting critical temperature T c 
and below a crossover temperature T* . The temperature 
T* decreases with increasing doping and merges eventu- 
ally to T c at (or slightly above) optimum dopingB. The 
pseudogap phase of underdoped cuprates is best char- 
acterized rk>Yi angle resolved photoemission spgictroscopy 
(ARPES)u □ and by tunneling experimentsErEl, which 
probe the single-particle excitation spectra directly. The 
pseudogap opening below T* corresponds to a suppres- 
sion of the low-frequency differential conductance (which 
is connected to the density of states) measured by tunnel- 
ing, and to a leading-edge shift of the spectral intensity 
(which is connected to the spectral function via the Fermi 
distribution and a dipole matrix element) measured by 
ARPES. As clearly shown by ARPES, the pseudogap is 
tied to the Fermi surface and its two-dimensional wave- 
vector dependence is strongly anisotropic, resembling a 
d x 2_ y 2 harmonic. Both ARPES and tunneling experi- 



ments suggest that the pseudogap evolves smoothly into 
the superconducting gap as the temperature is lowered 
from T* to T c . 

The d-wave-like wave-vector dependence of the pseu- 
dogap, its continuous evolution into the superconducting 
gap below T c , and its tying to the Fermi surface suggest 
that the pseudogap phase could be a precursor of the 
superconducting phase (at least in _Bi-based compounds 
for which a detailed ARPES analysis of the pseudogap is 
available). According to this interpretation, the crossover 
temperature T* acquires the meaning of the tempera- 
ture at which fluctuating pairs start forming without co- 
herence, the latter being not yet established owing to 
large fluctuations of the superconducting order parame- 
ter. Upon lowering the temperature, the coherence be- 
tween pairs is eventually established and superconduc- 
tivity appears. The occurrence of large superconduct- 
ing pair fluctuations in cuprates is related to the quasi- 
bidimcnsionality induced by the layered structure, as well 
as to the short coherence length £ of the superconduct- 
ing pairs, of the order of few lattice spacing (typically, 
Co -10-201). 

Within this scheme, the phase diagram of cuprates is 
interpreted in terms of a crossover from Bose-Einstein 
(BE) condensation of preformed pajcs,to BCS supercon- 
ductivity, as the doping is variedLnij. Heavily under- 
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doped cuprates are accordingly considered as supercon- 
ductors in a strong-coupling (BE) regime with T* 3> T c ; 
optimally doped and overdoped cuprates are instead 
more conventional superconductors in an intermediate- 
or weak-coupling (BCS) regime with T* ~ T c . The 
evolution from strong- to weak-coupling superconductiv- 
ity as the doping is increased is further supported by 
low-temperature ARPES and tunneling measurements 
in £?i-based compounds of the maximum superconduct- 
ing gap A (i.e., the gap at the M points of the Bril- 
louin zone), whereby A decreases as the doping is in- 
creased, with A ~ 60 — 70 meV in underdoped cuprates 
and Ao ™ 20 — 30 meV in optimally and overdoped 
cupratesBQ. Moreover, in underdoped Bi-based cuprates 
Ao is larger than the band-width along the M—Y(X) di- 
rections, suggesting that at least near the M points states 
with bosonic character can be formed. Recent high- 
resolution ARPES measurements on La2-xSr x CuOi at 
T = 11K also indicate that the superconducting gap near 
the M points increases as the doping is decreased, and 
smoothly evolves, into the normal phase pseudogap at 
very low dopingJiS. The doping dependence of the gap is 
thus likely to be an universal feature of cuprate super- 
conductors. 

Recent high-resolution ARPES experiments in Bi2212 
further suggest that a crossover from weak to strong cou- 
pling can even_be found along the Fermi surface (FS) 
at fixed dopingtH. Fermionic states near the nodal (N) 
points of the FS (namely, the points along the T — Y(X) 
directions where the band crosses the Fermi level and the 
d-wave gap vanishes) appear to be weakly coupled, while 
states near the M points are strongly coupled and could 
display bosonic character. Moving from N toward M 
points along the FS, a continuous crossover from weakly 
to strongly coupled states should accordingly be ob- 
served. In agreement with the expectation that increas- 
ing the coupling should cause an increase of the width 
of the spectral peaks, explicit support to the wave-vector 
induced crossover along the FS is obtained, for instance, 
from Fig. 2 of Ref. |7| The ARPES spectral intensities for 
an optimally doped Bi2212 sample reported in that figure 
show, in fact, that the width of the quasi-particle peak 
in the normal phase (as a function of both wave vector 
and frequency) increases along the FS, as one moves from 
the N toward the M point. In particular, near the M 
points the frequency distribution of the spectral intensity 
is broad and flat without any observable peak, while a 
broad peak feature is present in the wave- vector distribu- 
tion. In addition, in Bi2212 at optimum doping the band 
dispersion near the M points along the M — Y(X) direc- 
tion is rather narrow (of the order of 50 meV) , while the 
band dispersion along the T — Y(X) directions is consid- 
erably larger (of the order of 400 meV). For all cuprates 
for which ARPES measurements are available, the Fermi 
velocity vf is also anisotropic along the Fermi surface, 
with vf(N)/vf(M) ~ 3.E3 As a consequence, fermionic 
states near the M points are locally associated with a 
small Fermi velocity and strong coupling (hot fermions); 



while fermionic states near the N points are locally as- 
sociated with a large Fermi velocity and weak coupling 
(cold fermions). To account explicitly for the different 
properties about the M aneLJV points, a two-gap model 
has been recently proposed^. 

In the present paper, we investigate the evolution of 
the spectral function from the weak- to strong- coupling 
regimes in a systematic way, to compare with the evo- 
lution of the spectral function in cuprates by varying 
doping and wave vector. More specifically, we aim to ac- 
count for the character of the fermionic states near the M 
points (where bosonic states can be formed upon reduc- 
ing the doping due to the hot character of these states) 
and to follow the wave- vector induced crossover along the 
Fermi surface. The local character of the fermionic states 
in wave-vector space further enables us to use a simple 
isotropic attraction between electrons, which gives rise in 
the superconducting state to a gap with s-wave symme- 
try. 

Two different (albeit related) kinds of approaches for 
the pseudogap state can be identified within the pair- 
ing scenario. On the one hand, owing to the short co- 
herence length and the large. yaJpe. of the superconduct- 
ing gap about the M pointsoii3oE3, the superconduct- 
ing phase of underdoped cuprates is interpreted as in- 
termediate between a BCS state with extended pairs 
and a Bose-Einstein condensate with preformed (local) 
pairs. Within this view, due to strong- or intermediate- 
coupling effects, pairing correlations survive well above 
T c and determine a pseudogap opening when coupled to 
the fermions. On the other hand, the second approach 
emphasizes the relevance of phase fluctuations of the su- 
perconducting order parameter, owing to the low value of 
the plasma fxeouency and the quasi-bidimensionality of 
the cupratescJ'Ell. Within this view, the amplitude of the 
local order parameter is established at T* , even though 
phase coherence and hence long-range superconductivity 
occurs at the lower temperature T c . 

The approach we follow in this paper belongs to the 
first group of the pairing scenario. Specifically, we inves- 
tigate in a systematic way the role played by pair fluc- 
tuations in the pseudogap opening, following the BCS 
to Bose-Einstein crossover from weak to strong coupling. 
To this end, we introduce a simplified microscopic model 
representing a 3D continuum of fermions mutually inter- 
acting via an attractive contact potential, which can be 
parametrized in terms of the scattering length. This 3D 
model allows us to considerably simplify the numerical 
calculations as well to obtain analytic results (at least in 
some limits), yet preserving the qualitative features ob- 
tained for more realistic modela-like the two-dimensional 
negative-U Hubbard modelE3t3. 

We examine initially the two-particle propagator in 
the particle-particle channel, and evaluate the pair- 
fluctuation propagator T(k, uj) as a function of wave vec- 
tor k and frequency u>. We further analyze the single- 
particle propagator, and evaluate the self-energy E(k, cj) 
and the spectral function A(k, uS) within the non-self- 
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consistent t-matrix approximation. In the strong- and 
weak-coupling regimes, we discuss analytic forms for the 
self-energy, and comment on the main differences in the 
line shape of A(k, u>) between the two regimes. The spec- 
tral weight of the incoherent peak that appears in A(k, uj) 
and the temperature dependence of the chemical poten- 
tial are also discussed. In the intermediate (crossover) 
region (where analytic calculations are not feasible) only 
numerical results are presented. Our findings of different 
characteristic features occurring for A(k, oj) in different 
coupling regimes are then organized in a systematic way, 
and a criterion to distinguish by ARPES experiments 
whether an interacting fermion system is in the strong-, 
intermediate-, or weak-coupling regime is discussed. In 
the strong-coupling regime, we find it is appropriate to in- 
troduce two different crossover temperatures (T* and T *) 
to describe the peculiar evolution of the spectral func- 
tion for increasing temperature. We also show how these 
two temperatures merge to a single crossover tempera- 
ture (T*) as the coupling is decreased. A detailed com- 
parison of ARPES experiments with our systematic anal- 
ysis of the spectral function in different coupling regimes 
is eventually attempted. Although it might at first ap- 
pear that our model could not be directly applied for 
comparison with ARPES experiments in cuprates, this 
comparison is attempted by invoking the wave-vector- 
induced crossover mentioned above. Finally, we provide 
a heuristic method to rationalize the gross feautures of 
the evolution of the spectral function in terms of the BCS 
spectral function, whereby the BCS gap at T = is re- 
placed by the value of the pseudogap at (or above) T c and 
the Fermi energy is replaced by the renormalized chem- 
ical potential. Even though some parts of our analysis 
and results have been already presented (albeit for dif- 
feren t models .and/or with different methods) in previous 
worifflM, our approach should be regarded as more 
systematic and complete than others. 

The plan of the paper is as follows. In Section II we in- 
troduce the microscopic model and discuss the relevant 
equations for the spectral function and related quanti- 
ties. In Sections III and IV we report the results for two- 
and single-particle properties, respectively, discussing in 
a systematic way the evolution of the pair-fluctuation 
propagator and of the spectral function from the super- 
conducting fluctuations regime (weak coupling) to the 
bosonic limit (strong coupling). In Section V we present 
a detailed comparison of our results with ARPES exper- 
iments. Section VI gives our conclusions. 



II. RELEVANT EQUATIONS FOR THE 
SPECTRAL FUNCTION AND RELATED 
QUANTITIES 

In this Section, we set up the relevant equations to fol- 
low the evolution of the single-particle spectral function 
and the two-particle fluctuation propagator from weak 
to strong coupling. To this end, we consider a system 



of fermions embedded in a three-dimensional continuum 
and mutually interacting via an effective short-range at- 
tractive potential vo5(r — r') of strenght vq, where vo is a 
negative constant. For the 3D continuum model we are 
allowed to take the limit of a strictly short-range interac- 
tion (~ 6(r— r')), thus relating to the fermionic scattering 
length of by a suitable regularization procedure. Knowl- 
edge of the detailed form of the fermionic interaction is, 
in fact, not required for studying the main features of the 
evolution from weak to strong coupling. The many-body 
diagrammatic structure for the single- and two-particle 
Green's functions gets in this way considerably simpli- 
fied, while preserving the physical effects of pseudogap 
openingca. (For a detailed discussion of this model, see 
Ref . |f) 

For the attractive fermionic interaction of interest, the 
scattering length of changes from being negative (when 
the two-body problem fails to support a bound state) to 
being positive (when the bound state is eventually sup- 
ported by increasing the interaction strength), and di- 
verges when the coupling strength suffices for the bound 
state to appear. The dimensionless parameter kpciF 
(where kp is the Fermi wave vector) thus locates the side 
of the crossover one is examining and how close to the 
crossover region one is. Specifically, kpap is small and 
negative in the weak-coupling regime, diverges in the in- 
termediate (crossover) regime, and eventually becomes 
small and positive in the strong-coupling regime. For 
this reason, driving the crossover by varying kp while 
keeping ap fixed requires one to change discontinuosly 
the sign of ap at the value {kpap)^ 1 — 0. 

The diagrammatic scheme we consider is based on the 
non-self-consistent t-matrix approximation, constructed 
with "bare" single-particle Green's functions (with the 
inclusion, however, of the dressed chemical potential and 
of an additional constant energy shift (to be discussed 
below) which is relevant to the symmetry of the spec- 
tral function). This choice embodies the physics of the 
pseudogap state, because in the weak-coupling regime it 
describes the Ginzburg-Landau superconducting fluctu- 
ations above T]Ej, while in the strong-coupling regime 
it describes the formatkm of noninteracting bosons 
(fermionic bound states )Ea. 

The set of relevant equations for the two-particle 
Green's function in the particle-particle channel (i.e., the 
pair- fluctuation propagator), the single-particle Green's 
function, and the self-energy is the following: 

S (k,^ = -r£j|^r(°)( q ,a,) 
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G-\k,uj n ) = G(°r'(k, Wn ) - (53(k,w„) - S ), (3) 
,= 2T£e* + -/|^G(k,u,„). (4) 

Here, G^°^(k, w ra ) is the "bare" fermion propagator given 



A(k,w) = 



— ImE(k, u)/ir 



by G(°)"(k, W „) 



£(k) (£(k) = k 2 /(2m) - (j,' 



being the free-particle dispersion measured with respect 
to the renormalizcd chemical potential fj/ = [i — So, 
where fi is the physical chemical potential and E the 
constant self-energy shift mentioned above), m is the 
free-fermion mass, and u) n = irT(2n+ 1) (n integer) and 
f2„ = 2ixTv (y integer) are, respectively, fermionic and 
bosonic Matsubara frequencies at temperature T. The 
chemical potential is eliminated in favor of the density n 
via Eq. (4). The constant self-energy shift So is given 
by ReS(fc = fe p /,w = 0,T ~ T*) where = 
It turns out that this shift is non-negligible only in the 
weak- to intermediate-coupling regime, where it is al- 
most temperature independent from T to T*. [In this 
coupling regime, T* = T * = T£ is the temperature at 
which the pseudogap disappears.] In practice, we will 



take Sq = ReS(fc = fe p 



0, T = T c ) in the weak- 



to intermediate-coupling regime, while we shall neglect 
So altogether in the intermediate- to strong-coupling 
regime. Inclusion of this self-energy shift amounts to 
a partial self-consistency dressing of the single-particle 
Green's functions G^°\ and corresponds to a complete 
description of the high-temperature region. 

After analytic continuation to the real frequency 
axisEj, the imaginary part of the retarded self-energy can 
be written at this order of approximation in the form: 



ImS(k, uj) 



d 3 q 

Imr(°)(q,w + C(q-k)) 



[&(" + f(q-k)) + /(e(q-k))] 

(5) 



where f(x) = l/(e /3x + 1) is the Fermi distribution and 
b{x) = l/{e Px -l) is the Bose distribution, with f3 = 1/T. 
Hereafter, analytic continuations are meant to produce 
retarded (R) functions. Once the imaginary part of the 
self-energy is evaluated as above, its real part ReS(k, uj) 
is obtained via a Kramers-Kronig transform. The real- 
frequency formulation (5) allows for high accuracy of 
the numerical calculations, and avoids the problems of 
dealing numerically with analytic continuation from the 
imaginary frequency axis. 

The spectral function A(k, u>) for the single-particle 
fermionic excitations of interest is obtained from the 
imaginary part of the retarded Green's function G R (k, uj) 
via the relation 



A(k,w) = -~ImG R (k,uj 



(6) 



In terms of the real and imaginary parts of the self- 
energy, the spectral function A(k, uj) has the form: 



e(k)-ReS(k,^)+S ) 



(ImS(k,w)) 2 
(7) 



The set of equations (l)-(5), together with the defi- 
nition of the spectral function (7) and the prescription 
for analytic continuation to real frequencies, allows us to 
study in a systematic way with limited computational ef- 
fort two- and single-particle properties over a wide range 
of parameters (namely, coupling, density, and tempera- 
ture), following the crossover from weak to strong cou- 
pling. 

A few additional comments are in order at this point 
about the choice (1) of the self-energy. In the weak- 
coupling regime and at high enough temperature, the ex- 
pression (1) represents the leading term of a low-density 
expansion foe, a Fermi system even when the interaction 
is attractiveEZl. Upon approaching T c , the expression (1) 
can alternatively be interpreted as representing the cou- 
pling of a bare fermion with pairing fluctuations. In the 
strong-coupling regime, the expression (1) represents in- 
stead a "free" boson coupled to a bare fermion, and is 
known to produce a shadow-band structure in the spec- 
tral function at negative frequencies. We note that, to be 
consistent with a low-density approach, the "free" boson 
in Eq.(l) should be dressed with suitable self-energy cor- 
rections for composite bosons (as discussed in Ref. p7| ), 
at least in the not-too-extreme strong-coupling regime 
where the residual interaction between the composite 
bosons remains active. We defer this study to a future 
work, and regard the present approach with a "free" bo- 
son only as a preliminary but essential step to establish 
the qualitative behavior of the spectral function in a sys- 
tematic way, for all coupling regimes and over a wide 
temperature range. 

Finally, we recall that inclusion of full self-consistency 
in the single-particle Green's functions entering Eq.(l) 
can be safely dismissed (at least in the weak- and strong- 
coupling regimes, where the use of expression (1) can 
be justified to start with), as discussed in Ref. 27 In 
contrast to the non-self-consistent t-matrix approxima- 
tion used in Eqs. (2-4), the self-consistent t-matrix 
uses full self-consistent Green's functions but does not 
include vertex corrections in the self-energy. As re- 
marked in Ref. [2^, the different levels of approxima- 
tion for vertexes and single-particle Green's functions 
may then lead to unphysical results for the pseudogap 
in the spectral function, in a similar way to, what hap- 
pens for the 2D repulsive Hubbard mode£3. The im- 
portant point to be emphasized is that, provided G*- -* in 
Eqs.(l) and (2) contains the dressed chemical potential 
obtained from Eqs. (3) and (4), this set of equations in- 
terpolates smoothly between weak- and strong-coupling 
limits and provides a reasonable description of both lim- 
its. Moreover, we emphasize again that our interest in 
the (three-dimensional) continuum model originates es- 
sentially from obtaining analytic results in the continuum 
case, from which the main features of the spectral func- 
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tion can be readily extracted. A comparison between our 
systematic results for the continuum model and those 
available for the lattice model will be made in Section 
VI, with the outcome that the main qualitative features 
remain the same in the two models. 



III. PAIRING FLUCTUATIONS FROM WEAK 
TO STRONG COUPLING 

In this Section, a systematic study of the crossover from 
the superconducting fluctuation (weak-coupling) regime 
to the bosonic (strong-coupling) regime for the pair- 
fluctuation propagator is reported via both analytical 
and numerical calculations, for a wide temperature range 
above T c . This study is preliminary to the discussion of 
the spectral function via the self-energy (1), presented in 
the next Section. 

The pair-fluctuation propagator r(°)(q, Sl„) within the 
non-self-consistent t-matrix approximation is given by 
Eq.(2). We shall examine, in particular, its wave- vector 
and frequency dependence for all coupling regimes. 

Being interested in normal-phase properties of the 
fermionic system, knowledge of the superconducting crit- 
ical temperature T c is required at the outset to insure 
that T > T c . To identify the critical temperature T c , 
we rely on the condition that the fluctuation propagator 
has a pole at T c for vanishing wave vector and frequency, 
namely, r( ) _1 (q = 0, tt v = 0;/x,T = T c ) = 0. This 
condition (known as the Thouless criterion for the super- 
conducting instability) taken alone is equivalent to the 
BCS equation for the critical temperature in the weak- 
coupling limit. In addition, when coupled to the density 
equation (4) to fix the chemical potential, it yields the 
value of the Bose-Einstein condensation temperature in 
the strong-coupling limit. 

The Thouless criterion provides a temperature- 
dependent critical value for the chemical potential 
fJ-c(T c ) = ^(n,T c ), with /i(n,T) obtained from the den- 
sity equation. (Recall once more that, for a proper 
description of the crossover from weak to strong cou- 
pling, it is essential to let the chemical potential ad- 
just itself with coupling at given density.) In Fig.l 
the critical chemical potential and the normal-phase 
chemical potential are reported for different values of 
(fcFdf) -1 in the strong- to intermediate-coupling regime. 
We have verified that in the high-temperature limit 
the chemical potential /i(n,T) tends to its classical 
value fi(n,T) = 1.5Tln(1.898 T B e/T) for all densities 
(Tbe = 3-3ln 2 J 3 /rux being the Bose-Einstein condensa- 
tion temperature). E^l 
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FIG. 1. Critical chemical potential ii c (T c ) (full line) and 
normal-phase chemical potential fi(n, T) for different values 
of the density expressed in terms of (^fcsf) - (broken lines), 
as functions of temperature. The values of (A^af) -1 are 
1.33,0.96,0.76,0.60, and 0.52 from bottom to top. 

Note that in the strong-coupling regime, i.e., at low 
density and for (3\[i\ ^> 1, [i approaches the value —eq/2, 
where eo = (maf*) is the binding energy of the asso- 
ciated two-body problem. In the intermediate-coupling 
regime, /i initially increases as the temperature is in- 
creased, reaches a maximum, and eventually decreases, 
tending to the classical behavior. The temperature where 
the maximum is located turns out to be smaller than 
the binding energy eo, and does not relate to specific 
changes of single-particle properties in the strong- or 
intermediate-coupling regimes. In the weak-coupling 
regime, on the other hand, this temperature turns out to 
coincide with the crossover temperature T* where pair 
fluctuations become manifest and a pseudogap opens. 

In the weak-coupling regime, the presence of the max- 
imum in [i(T) is connected with the fermionic degrees of 
freedom, and in particular with the opening of a pseu- 
dogap at the Fermi surface which tends to depress the 
chemical potential upon lowering the temperature (in a 
similar fashion to what happens when a real gap opens 
in a BCS superconductor below T c ). The presence of a 
maximum of /i(T) (with the ensuing non monotonic be- 
havior of /i(T)) is clearly observed in Monte Carlo sim- 
ulations of the 2D attractive (s-wave) Hubbard model 
in the intermediate-coupling regime (cf. Fig. 6 of Ref. 

), while the presence of the maximum is somewhat de- 
batejcLwhen the self-consistent t-matrix approximation is 
usedEaa. In the strong-coupling regime, the fermionic 
degrees of freedom are exponentially suppressed accord- 
ing to /(£) ~ exp(— and the above maximum is 
progressively shifted toward zero temperature for increas- 
ing j3\n\, thus recovering in the extreme strong-coupling 
limit the behavipx-oL a free Bose gas via the relation 
2/i = — eo + /isoE3i3. We have verified that, in the 
strong- and intermediate-coupling regimes, the tempera- 
ture at which n(T) reaches its maximum does not relate 
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with the temperature at which the pseudogap opens. 

The equation for the density, together with the con- 
dition for the critical chemical potential (i(T) = n c (T), 
yields the value of the critical temperature T c for the 
superconducting instability. The critical temperature T c 
and the BCS mean-field critical temperature Tbcs (the 
latter obtained formally from the same equations defin- 
ing T c but with the bare single-particle Green's function 
q(o) replacing G in Eq.(4)) are reported in Fig. 2 as func- 
tions of the parameter (kpaF)~ l (recall that in the weak- 
coupling limit Tbcs — l-67e.p exp(7r/2fc^af ) (ap < 0)). 
(Both temperatures have been conveniently normalized 
to the Bose-Einstein condensation temperature Tb e eval- 
uated at the same density. Note that in the 3D contin- 
uum model Tbe is of the same order of e^.) The re- 
sults of Fig. 2 can be compared with the calculation of 
T c within the non-self-consistent t-matrix approximation 
reported in Ref . We mention that the inclusion of the 
self-energy shift So in the weak- to intermediate-coupling 
regime adopted in the present paper, slightly increases 
the value of the critical temperature with respect to the 
results of Ref. |35]. 
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FIG. 2. Critical temperature T c (full line) and BCS 
mean-field critical temperature Tbcs (broken line) as func- 
tions of (^fof)" 1 (both temperatures are normalized to the 
Bose-Einstein condensation temperature Tbe evaluated at 
the same density). 

We have verified that, in the strong- and intermediate- 
coupling regimes, the mean-field temperature Tbcs 
about coincides with the temperature at which the bare 
(fermionic) contribution uq to the total density n equals 
the (bosonic) contribution 5n due to interaction effects, 
i.e., n (T = T B cs) - Sn(T = T BC s) - n/2. This result 
permits us to identify Tbcs as the crossover temperature 
where preformed pairs start to form. [The connection 
between Tbcs and the characteristic crossover temper- 
ature^) of the spectral function will be made in Sec- 
tion IV.] Note from Fig.2 that for (kpap)- 1 < -1 the 
critical temperature approaches the BCS mean-field crit- 
ical temperature, indicating that the fermionic system 



is in the weak-coupling regime. For (fcpai?) -1 > I the 
critical temperature is instead close to the Bose-Einstein 
condensation temperature, indicating that in the strong- 
coupling regime the fermionic system is equivalent to a 
system of non-interacting bosons. The strong-coupling 
limit is thus effectively reached for not too large values 
of the parameter (kpap) -1 S3 

In Fig. 3 we report for convenience the relation between 
kF^pair and {kpap) -1 as obtained from the analytic so- 
lution of Ref. [58|, £ P air being the average pair size at 
T = 0. This plot is especially useful to compare our re- 
sults for the two- and single-particle properties (which are 
expressed in terms of (kpap)" 1 ) with the phenomenology 
of cuprates, for which some estimates of the parameter 
kp^pair in different doping regimes are available!! 
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FIG. 3. Phenomenological parameter fc^Cpair as a function 
of the coupling parameter (fcaF) -1 , from weak to strong 
coupling. 

Specifically, at optimum doping the parameter kF^pair 
takes roughly values between 6 and 10, and its value de- 
creases for decreasing doping (mainly because the super- 
conducting gap at T — increases and the Fermi energy 
decreases approaching the insulating phase). We may 
reasonably consider kp^pair ~ 1 as a lower-bound value 
for underdoped cuprates, especially if we consider it as 
a local quantity about the M points. Accordingly, the 
coupling parameter (fc^a^) -1 in the optimum and un- 
derdoped regimes for cuprates lies approximately in the 
range —1.7 < (kpap)^ 1 ^ 0.5, as indicated in Fig. 3. 

Having determined the thermodynamic quantities 
/Lt(n, T) and T c , we pass now to calculate the pair- 
fluctuation propagator (2). From a physical point of 
view, pairing fluctuations have essentially different char- 
acter in the strong- and weak-coupling regimes. While to 
evaluate the self-energy numerically from Eq.(5) knowl- 
edge of r(°) is required over a wide range of wave vectors 
and frequencies, to characterize the evolution of T^ ' from 
weak to strong coupling it is sufficient to consider the ex- 
pansion of its inverse in poweps, pf the wave vector q and 
the Matsubara frequency OyJld'ta 



G 



r (orl (q,^) 



6|q| 2 + din,, 



(8) 



with d — d\ + id2Sgn(p, v ). Here, (a, 6, di, c^) are real 
coefficients which are coupling, density, and temperature 
dependent. 

In Fig. 4 we report the ratio of the imaginary (da) to 
the real (d\) part of the frequency coefficient in Eq.(8) 
at T c as a function of (fcfaf) -1 . 



CM 




(k F a F )- n 

FIG. 4. Ratio between the imaginary (d2) and real (di) 
part of the frequency coefficient of the inverse pair-fluctuation 
propagator at T c as a function of (^faf) -1 - 

In the strong-coupling limit, d\ ~ — m 2 ap /(8ir) and 
d% ~ 0, with the pair- fluctuation propagator acquir- 
ing the polar structure of a bosonic Green's function. 
In the weak-coupling limit, on the other hand, d\ ~ 
-(T C /E F ) 2 < I and d 2 = -N tt/(8T c ), where iV is the 
density of states (per spin component) at the Fermi en- 
ergy Ep, with the pair-fluctuation propagator acquiring 
the diffusive Ginzburg-Landau structures 
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FIG. 5. Coefficient b of the |q| 2 term in the inverse 
pair-fluctuation propagator at T — T c as a function of 

(fcflF)" 1 . 

When |di| > |d 2 |, fluctuating pairs become propagat- 



ing (albeit with a damping), and eventually acquire the 
(bosonic) character of undamped preformed pairs. 

In Fig. 5 the coefficient b of the |q| 2 term in Eq.(8) is 
plotted as a function of (kpap)~ l . 

In the strong-coupling regime (i.e., for (kpap) _1 > 1), 
b coincides with its strong-coupling value map / (32tt), 
where 2m is the mass of the composite boson. In the 
weak-coupling regime (i.e., for l/kpap < — 1), b = 
N 7 Q(3)v 2 F / {4&ir 2 T c 2 ) is proportional to the square of the 
zero-temperature correlation length £o (where vf is the 
Fermi velocity and £(3) ~ 1.202 is the Riemann zeta 
function of argument 3), and tends to diverge in the ex- 
treme weak-coupling limit. 

Finally, the coefficient a in Eq.(8) provides a mass to 
the pair-fluctuation propagator. In the weak-coupling 
limit a = No ln(T/T c ) vanishes at T c , while in the strong- 
coupling limit we can write a = —m 2 ap[j,B(T)/ (8ir) 
where /xb(T) is the bosonic chemical potential, which 
we pass now to discuss. 

According to the above analysis, we have verified that 
in the strong- coupling regime (when the conditions /i < 
and f3\n\ ^> 1 are satisfied) the pair- fluctuation propaga- 
tor evaluated numerically acquires the polar structure of 
a free-boson Green's functionL^fEZI: 



r (0) (q,^) 



87T / \m 2 ap) 



ifi„ - q 2 /(4m) + fi B 



(9) 



with mass 2m, a quadratic dispersion as a function 
of wave vector, and bosonic chemical potential /j,b(T), 
which reduces to 2p,(T) + e m the extreme strong- 
coupling limit. This implies that, for real frequencies, 
the imaginary part of r*- -* (which enters the calculation 
of the imaginary part of the self-energy via Eq.(5)) is 
proportional to a delta function in the strong-coupling 
limit: 



Imr (0) (q,w) 



8tt 2 
m 2 ap 



S(LU-q 2 /(4m) + f i B ). (10) 



We have verified that Eq.(10) remains approximately 
valid in the intermediate- coupling regime toward strong 
coupling (when (3\[i\ ~ 1 and fj, < 0), while only in the 
strong-coupling limit \x b {T) therein reduces to the chem- 
ical potential no(T) of an ideal Bose gas, with the char- 
acteristic temperature dependence 



M o(T) = -1.22 



(T - T BE ) 2 
Tbe 



(11) 



at low enough temperature. Specifically, we have ver- 
ified that a delta-function contribution to Imr(°) (q, uj) 
appears when (kpap)^ 1 > (corresponding to d 2 /di < 1 
in Fig. 4). This contribution, which is initially present for 
large values of q, extends progressively to smaller values 
of q for increasing coupling, reaching eventually q = 
when the chemical potential becomes negative. By fur- 
ther increasing the coupling, the delta-function contribu- 
tion to Imr(°)(q, u>) becomes increasingly prominent and 
the asymptotic expression (|l(]) is progressively reached. 
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This delta-function contribution, associated with the 
formation of a bound state with bosonic character, is 
responsible in the strong-coupling limit for the opening 
of a real gap in a broad temperature range above T c , as it 
will be shown in the next Section. Actually, in addition 
to the delta function, one also finds a finite contribution 
to ImT(°) for ui > 2 | fj, |. When inserted in Eq. (5), this 
contribution leads to an exponentially vanishing ImE for 
u)>\(jl\. 

Finally, in the weak-coupling regime, the pair- 
fluctuation propagator recovers the Ginzburg-Landau 
diffusive form. Near the critical temperature, its expres- 
sion for small wave vectors and frequencies is accordingly 
given by: 



r<°>(q,n„) = 



1 



iVo(£ + 7/|q| 2 +7|^|) 



(12) 



Here, e = ln(T/T c ) is the mass term of the propagator, 
V ~ "|/T C 2 represents the stiffness of the superconduct- 
ing fluctuations with a proportionality coefficient which 
depends on dimensionality (7£(3)/(487r 2 ) being its value 
in 3D), while 7 = 7r/(8T c ) is related to the lifetime of 
the fluctuating pairs (which do not obey Bose statistics). 
In this limit, Imr(°'(q = 0, ui) diverges like 1/ui only at 
the critical temperature T = T c ; as a consequence, the 
pseudogap region induced by the diffusive pair fluctua- 
tions will be present only in a rather narrow tempera- 
ture range. No delta function contributes in this regime, 
but Imr(°)(q = 0, ui) has a broadened peak structure for 
small enough q. 

To summarize, the main effect of increasing coupling in 
ImrW (q, ui) is the appearance of a peak structure (delta 
function) at finite frequencies, whose area grows with 
(fcpflF)" 1 - In the strong-coupling regime, for q = and 
T > T c the real part of r(°)(q, w)^ 1 vanishes at a fi- 
nite frequency, corresponding to a pair resonance. This 
resonance disperses as q 2 . In the weak-coupling regime, 
the real part of r'°^(q, ui)~ x is small only in the critical 
region, and vanishes only at the critical temperature. In- 
creasing the coupling from weak to strong, the frequency 
dependence of ImT^°^(q = 0,u>) evolves from being an- 
tisymmetric with respect to w = to an asymmetric 
structure. This evolution confirms previous results for 
r(°)(q = 0,w) reported in Ref. 

In the next Section, we will show how the peak struc- 
tures of r(°)(q, ui) affect the single-particle self-energy 
and hence the spectral function, giving rise to a (pro- 
nounced) suppression of the low-energy spectral weight, 
namely, to a pseudogap. 



IV. SPECTRAL FUNCTION FROM WEAK TO 
STRONG COUPLING 

In this Section, we study the single-particle excitations 
for fermions coupled to pair fluctuations above the criti- 
cal temperature. The spectral function A(k, ui), obtained 



by solving the set of equations (l)-(7), is analyzed in a 
systematic way as a function of coupling and tempera- 
ture, thus following its evolution from weak to strong cou- 
pling. In this way, characteristic features of the spectral 
function as a function of frequency and temperature will 
be evidenced in all coupling regimes. We shall analyze 
separately the cases when the chemical potential lies be- 
low the bottom of the single-particle band (fi < 0) corre- 
sponding to the strong- to intermediate-coupling regime, 
and when the chemical potential lies inside the single- 
particle band (zx > 0) corresponding to the intermediate- 
to weak-coupling regime. 



A. Strong- to intermediate-coupling regime 

In the previous Section, we have verified that a 
delta function appears in ImF^ (q, ui) starting from the 
intermediate-coupling regime when (kpap) 1 > (for 
temperatures such that f3\/i\ > 1). In particular, we have 
verified that in the strong- coupling limit (where /1 < 
and P\fx\ 3> 1) the pair- fluctuation propagator coincides 
with a free-boson Green's function with mass 2m (cf. 
Eq.(9)). In this limit, the imaginary part of the pair- 
fluctuation propagator reduces to a delta function and 
the self-energy can be evaluated analitically. Inserting 
Eq.(10) into the general expression (5) for the imaginary 
part of the self-energy, the following form results: 



Im£(k = 0,w) 



-2(4m) 3 / 2 y/uj th - ujQ(uj th - ui) 



m^aF 



O0{uith-U+\^B |) _ I 



(13) 



where uith = /i — /is is a threshold frequency and Q is 
the unit step function. [In the strong- to intermediate- 
coupling regime, when the chemical potential is below the 
bottom of the free-fermion band (// < 0) and k 2 / (2m) <C 
eo, the self-energy and hence the spectral function are 
almost independent of wave vector. In this case, k = 
can be taken as a representative value, as we did in 
Eq.(13).] 

Note that the frequency dependence of ImS is strongly 
asymmetric about its minimum at u) ~ uith — |Mb|. Note 
also that ImS (and hence ReS obtained via Kramers- 
Kronig transform) has a nontrivial temperature and 
frequency dependence, showing strong deviations from 
Fermi-liquid behavior. In the regime where /3|/Ab| << 1 
(i.e., T ~ T c ), three different behaviors of lm£(0, ui) can 
be specifically identified on the frequency axis: (i) For 
uith - I Ms I < ui < u) th , ImS(0,w) ~ -^Uth ~ ui/P\/j, b \; 
(ii) For uth - P < u) < uith - |a*b|, ImS(0,w) — 

— (fly/uJth - u;)^ 1 ; (Hi) For ui < uith - (3, lm£(0,w) ~ 

— exp(— P(u>th — w ))- Note that, in the strong-coupling 
limit, the imaginary part of the self-energy has a square- 
root divergence at u> = oj t h for T = T C . 

We have further verified numerically that strong de- 
viations from Fermi-liquid behavior are present in the 
strong- to intermediate-coupling regime in a wide tem- 
perature range above T c (while in the weak-coupling 
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regime non-Fermi-liquid behavior is found only in a nar- 
row temperature range above T c , as discussed in the next 
subsection). 

The above characteristic features of the analytic ex- 
pression (13) can be clearly identified in the numerical 
results for ImS(k = 0,ui) reported in Fig. 6 at different 
temperatures (for a specific coupling). The associated 
real part is shown in Fig. 7, where the straight lines to + fi 
are also reported for the same temperatures (and cou- 
pling) , with increasing temperature from top to bottom. 
At any temperature, the intersection of a given straight 
line with Re£(k = 0, u>) locates the position of the quasi- 
particle peak at to > 0. 
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FIG. 6. Imaginary part of the self-energy at k = as a 
function of frequency (in units of eo) at different temperatures. 
In this case, {kpapy 1 = 0.77 and T c /Tbe = 1.14. (Strong- 
to intermediate-coupling regime.) 
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FIG. 7. Real part of the self-energy at k = as a function 
of frequency (in units of eo) at different temperatures. In 
this case, (fcaf) -1 = 0.77 and T c /Tbe = 1.14. (Strong- to 
intermediate-coupling regime.) 

Note that for temperatures close to T c , three intersec- 
tions occur, with the most left intersection giving rise 



to the incoherent peak in A(k, a;) at negative frequen- 
cies (while to the central intersection there corresponds 
a strong suppression of A(k, uj)). At high enough tem- 
peratures, on the other hand, only a single intersection 
occurs. The corresponding spectral function for the same 
coupling and temperatures above T c is reported in FigJ|. 
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FIG. 8. Spectral function at k = as a function of fre- 
quency (in units of eo) at different temperatures. In this case, 
(kpap)' 1 = 0.77 and T c /T B e = 1.14. (Strong- to intermedi- 
ate-coupling regime.) 

The resulting spectral function has a strongly asym- 
metric structure with two peaks: The one at positive fre- 
quencies is rather narrow, coherent like, and has a large 
spectral weight (namely, the area enclosed by the peak); 
The one at negative frequencies is instead broad and has 
a small spectral weight. When the chemical potential 
/i is below the bottom of the band (as in the strong- 
coupling limit), the peak located at negative frequencies 
represents the incoherent peak generated by the interac- 
tion of the fermions with strong pair fluctuations (recall 
that, in this limit, the pair-fluctuation propagator has 
the polar structure of a bosonic Green's function). This 
incoherent peak is itself asymmetric, it becomes broader 
for increasing temperature, its spectral weight is density 
and coupling dependent (decreasing like (fci?ai?) 3 ), and 
its position depends mainly on the value of the chemical 
potential (which in turn depends on temperature). By 
increasing temperature, the chemical potential becomes 
progressively more negative (cf. Figjl]) and the peak po- 
sition ofA(k. — 0,uj) shifts accordingly toward positive 
frequencies. The broadening of the incoherent peak be- 
comes pronounced when the temperature is of the order 
of the binding energy (see, e.g., the case with T/T c = 3 
in Fig. 8). Note also that, for increasing temperature, the 
two peaks in A(k = 0, uj) get broadened in an asymmet- 
ric way (in contrast to the weak-coupling regime (see be- 
low), approaching which the broadening of the two peaks 
becomes progressively more symmetric). We have also 
verified that, in the extreme bosonic limit, the spectral 
function has the structure of two delta-like peaks sym- 







metrically located with respect to ui — (albeit with 
quite different spectral weights), which is generated in 
an asymmetric way by the narrowing of the incoherent 
peak at negative frequencies as the product kpap be- 
comes smaller and smaller. 

In Fig.^| the spectral function at k = is plotted as a 
function of frequency for different values of the parameter 
(kpap) -1 at T = T c . Note that these curves have been 
expressed in units cp, instead of eo, to get a more evident 
evolution with coupling. 
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FIG. 9. Spectral function at k = as a function of fre- 
quency (in units of ep) at T = T c for different values of 
the coupling (kpap) -1 . (Strong- to intermediate-coupling 
regime.) 

The spectral function has two well-separated peaks, 
with a real gap opening at an energy of the order of the 
binding energy of the pairs. By increasing the coupling, 
the spectral weight inside the gap is progressively sup- 
pressed, until in the extreme strong-coupling limit the 
step function in the imaginary part of the self-energy 
(13) makes the spectral weight to vanish identically in 
the range -\fi\ + \^b\ < u < ImI- 

Since photoemission experiments measure the intensity 
of photo-emitted electrons (that is, the spectral weight at 
negative frequencies) , no signal would be detected if both 
the incoherent and coherent peaks had moved to positive 
frequencies for increasing temperature. In this context 
(and in analogy with what is empirically done when in- 
terpreting photoemission measurements), it is natural to 
introduce a crossover temperature Tq at which the maxi- 
mum of the lower peak crosses zero frequency. Our anal- 
ysis shows, however, that at T * the spectral function still 
maintains a two-peak structure (cf. Fig.||), reflecting the 
sizable effects of the interaction between fermions and 
pair fluctuations. We are accordingly led to introduce 
a second crossover temperature T* > Tq, at which the 
upper and lower peaks of the spectral function merge 
just in one peak (in the sense that the incoherent peak 
is progressively absorbed by the coherent (quasi-particle) 
peak, even though the separation between the two peaks 



remains almost constant). [Photoemission measurements 
alone, however, would not be able to identify this second 
crossover temperature T* since the merging of the two 
peaks would occur at positive frequencies, a region which 
only inverse photoemission is able to probe.] 

In Fig. 10 the two crossover temperatures T * and T* 
(as obtained numerically from the above definitions) are 
reported as functions of the parameter (kpap) -1 , both 
temperatures being normalized with respect to the crit- 
ical temperature T c . The BCS mean-field critical tem- 
perature Trcs from Fig. 2 is also reported for compari- 
son. In the strong-coupling limit (when (kpap) -1 > 1) 
T* ^> Tq, T* being a large energy scale which, accord- 
ing to Fig. 10, in the strong-coupling limit is much 
larger than the binding energy eo- The difference be- 
tween T* and T * is reduced by decreasing (kpap) -1 , 
but only in the intermediate-coupling regime (i.e., when 
(kpap) -1 < —0.1) the two crossover temperatures al- 
most coincide (T* ~ Tq). In the weak-coupling regime, 
only a single crossover temperature can be identified 
(T* = Tq), since in this regime the chemical potential 
is almost equal to the Fermi energy and the two peaks 
of the spectral function are symmetrically located about 
zero frequency. Note finally that Tq about coincides with 
Trcs which was previously identified via an independent 
procedure. [When T * ~ T* we shall indicate both tem- 
peratures simply as T*.] 
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FIG. 10. Crossover temperatures Tq (dashed line) and Tj* 
(full line), and BCS critical temperature Tbcs (dashed-dotted 
line) as functions of (kpap) -1 ; all temperatures are normal- 
ized to the critical temperature T c of Fig. 2. The value of 
(kFap)^ 1 where the chemical potential changes sign is indi- 
cated by an arrow. 

B. Intermediate- to weak-coupling regime 

In the intermediate to weak-coupling regime, it be- 
comes essential to take explicit account of the constant 
shift Eo introduced in Section II. This shift has been iden- 
tified with the value of the real part of the self-energy (1) 
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taken at the wave- vector Ay (such that = and 

where the pseudogap turns out to be minimal), at zero 
frequency (about which the relevant range of the pseudo- 
gap phenomena is centered), and at a temperature close 
to T* (where the system recovers a Fermi-liquid behav- 
ior). An exact selection of T is, in practice, not required 
since Eo turns out to depend rather weakly on T in the 
intermediate- to weak-coupling regime. In this sense, we 
interpret Eo as a kind of Hartree shift, even though for 
our choice of the potential the true Hartree shift vanishes 
identically. We have consistently evaluated the constant 
shift Eo at |k| = fc^/,w = 0, and at the same temperature 
where the self-energy E(k, to) of Eq.(l) is calculated!^. 

The inclusion of the above constant shift Eo stems 
from the need of improving the single-particle Green's 
functions entering the construction of the self-energy (1) 
when approaching T c , only close to which pseudogap phe- 
nomena become appreciable in the intermediate- to weak- 
coupling regime. The choice of the self-energy (1) takes, 
in fact, into account fluctuation corrections only at the 
lowest order, a procedure which is certainly not com- 
pletely satisfactory when approaching the critical tem- 
perature where all sort of fluctuations corrections be- 
come important. To approach T c , one may try to improve 
the self-energy (1) by dressing the single-particle Green's 
functions therein with a constant self-energy insertion ap- 
propriate to the non-critical (temperature) region. On 
the other hand, the inclusion of the full self-consistent 
Green's function (without vertex corrections, however) 
leads to an overall depression o f nsy udogap phenomena 
and is not theoretically justifiedcjO. From a pragmatic 
point of view, we have verified that in the intermediate- 
to weak-coupling regime the pseudogap would open at 
negative frequencies (and not at w = 0, as expected 
from a simple physical intuition), if the constant self- 
energy shift Eo were not properly included. The pseu- 
dogap opening at negative frequencies would, in turn, be 
in contrast with Monte Carlo results and experimental 
findings. 

The characteristic behavior of the imaginary and real 
parts of the self-energy at |k| = k^i are shown in Figs. 11 
and 12, respectively, at different temperatures (for a 
given coupling). Note that the convexity of the curves 
ImE(fcp',w) about u> = is inverted with respect to the 
Fermi-liquid behavior, implying strong deviations from 
Fermi-liquid behavior also at moderate values of the cou- 
pling (i.e., such that a bound-state in the two-body prob- 
lem is not yet present). We have verified, however, that 
the Fermi liquid behavior is consistently recovered when 
the coupling is progressively decreased. 

In the weak-coupling limit and for temperature close 
to T c , an analytic approximation for the imaginary 
part of the self-energy can be obtained by inserting 
into Eq.(5) the weak-coupling expression of the pair- 
fluctuation propagator given by Eq.(12). At zero fre- 
quency and at the Fermi wave vector, the imaginary part 
of the self-energy acquires then the following expression 
in the limit T — > T r : 



ImE(k F ,cj = 0) 



6^ (TV) Tc, 
7C(3) WJ 2 n 



T-T c 



(14) 



which diverges upon approaching T c with a slow loga- 
rithmic rate. An expression analogous to (14) is also 
obtained at finite frequency (such that | lo |^cC ep) and 
T = T c , with the replacement of ln((T - T c )/T c ) by 
ln(| lo | /lo c ), where w c C £f is a suitable cutoff fre- 
quency. 

To test the validity of the above analytic approxima- 
tions, we may consider, e.g., the case of Fig. 11 for T/T c — 
1.001 and obtain from Eq.(14) the value lmT,/e F ~ -0.92 
for lo = 0. This estimate is indeed in good agreement 
with the numerical result reported in Fig. 11 (cf. the full 
curve therein), for which ImE/ei? ~ — 1. A fine tuning of 
the temperature very close to T c is, however, necessary 
to get a sizable increase of | ImE | due to the logarith- 
mic divergence in Eq.(14). For instance, to double the 
above value a temperature (T — T c )/T c = 10~ 6 has to 
be reached. In the 3D model here considered, the di- 
vergence of ImE is therefore not numerically detectable 
for all practical purposes. In addition, to test the va- 
lidity of the counterpart of Eq.(14) extended to finite 
frequency as explained above, we may consider the case 
of Fig. 11 for T/T c = 1.001 and two different frequencies, 
say, LO\jtF — 0.075 and Lo 2 /eF = 0.037. In this case, 
we obtain from our analytic approximation the value 
(ImE(aji) — ImE(w2))/eF = 0.093, which is rather close 
to the numerical result 0.106 as obtained from Fig. 11. 

The analytic approximation (14) [as well as its counter- 
part at T — T c and finite lo] need to be compared with the 
analytic form of E(k, lo) obtained in the weak-coupling 
limit by Ref. ^ within the same non-self-consistent t- 
matrix approximation adopted in the present paper. Ac- 
cording to Ref. |{J the diffusive form (12) of the pair- 
fluctuation propagator would yield 



E(k,w) 



A 2 

pg 



lo + £(k) + «7 



(15) 



where A pg is a parameter that depends on a wave- vector 
cutoff and 7 oc (T — T c ). This expression evidently does 
not reduce to Eq.(14) for lo — and T — » T c , nor to the 
counterpart of Eq.(14) for T = T c and finite lo. A few 
comments to clarify the origin of these discrepancies are 
then in order. 

The expression (|l5| ) has been derived more recently 
in Ref. where it was also extensively used to fit 
ARPES data for Si-based cuprates. According to Ref. 
[44] , Eq.(|T5|) results by manipulating directly the ex- 
pression (Tj for the self-energy in Matsubara frequency, 
whereby the finite value ttT of the smallest (fermionic) 
Matsubara frequency is exploited to make approxima- 
tions on the q-dependence of the integrand. Analytic 
continuation to the real frequency axis is then performed 
on the approximate result, yielding eventually the ex- 
pression ( |T5| ) above. This procedure is, however, ques- 
tionable, insofar as the very variable to be analytically 
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continued is used to set restrictions on the approximate 
form of the function (in this case, the q-dependence of 
the integrand). In our procedure, on the other hand, an- 
alytic continuation is performed at the outset [cf. Eq.(5)] 
and the relevant (controlled) approximations to get the 
approximate result (14) are introduced only afterwards. 

Note, in addition, that at | k |= fep the expression 
( |l5| ) produces two peaks symmetrically located about 
to = 0. This expression cannot, therefore, be used to 
fit the curves of A(k )Ji > , to) for the coupling values we are 
considering [cf. Figs. 13 and 16 below], whereby the sym- 
metry of the two peaks is recovered only in the extreme 
weak-coupling limit. In the analysis reported in Ref. ^4|, 
on the other hand, the experimental data are artificially 
symmetrized and the expression (^]) (together with an 
additional scattering rate iT\) is used to fit the ARPES 
data. We shall propose below an alternative phenomeno- 
logical fit to the curves of A(k tl /,uj), which is suggested 
by our numerical calculations. 

In our numerical calculations we have found that, at 
low enough temperature, there are three intersections of 
the curves (Re£(fc 



to) — So) with the straight line to 
(not shown in Fig. 12), with the two outer intersections 
giving rise to the two peaks of A(k^ , to) (see Fig. 13) while 
the central intersection corresponds to a strong suppres- 
sion of A(k fl ' , to) owing to the associated large value of 
ImE (A;^/, cj). By increasing temperature, on the other 
hand, only one intersection remains (resulting in only 
one visible peak in A(k^,u>) - see Fig. 13). 
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function of frequency (in units of ep) at different temperatures 
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The associated spectral function at |k| = fc M / is re- 
ported in Fig. 13 for the same temperatures and coupling 
of Figs. 11 and 12. 
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FIG. 13. Spectral function at |k| = fc M / as a function of 
frequency (in units of e F ) at different temperatures. In this 
case, with (/cfOf) -1 = —0.45 (T c /eF = 0.23). (Intermediate- 
to weak-coupling regime.) 

The spectral function obtained in the intermediate- 
coupling regime shows a well-developed two-peak struc- 
ture near T c with a minimum at zero frequency; yet the 
spectral weight distribution remains slightly asymmetric 
about zero frequency, even when approaching the crit- 
ical temperature. At zero frequency the spectral func- 
tion has a sizeable finite value, indicating that no real 
gap opens at the Fermi surface. Note from Fig. 13 that, 
upon increasing the temperature, the pseudogap fills in 
and closes at the same time, with the two peaks of the 
spectral function merging in just one peak at a crossover 
temperature T* ~ T * (which in this particular case is 
between 1.05T C and 1.08T C ). It is thus apparent that a 
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breakdown of the normal-state Fermi liquid occurs well 
before the system is in the preformed-pair limit. From 
the two-peak structure of .A (Ay , oj) in the intermediate- 
to weak-coupling regime, a pseudogap A pg could be em- 
pirically defined either as half the frequency separation 
between the maxima of the peaks, or as the separation of 
the maximum of the lower peak (at negative frequencies) 
from zero frequency. These two definitions coincide in 
the weak-coupling limit but slightly differ in the strong- 
coupling limit (see also Table II below) . Throughout this 
paper we will adopt the second definition, which is the 
most relevant for comparison with photocmission exper- 
iments, accessing only negative frequencies. 

In the intermediate-coupling regime, when the chemi- 
cal potential lies inside the fermion band and the Fermi 
surface is well defined, the wave-vector dependence of 
the spectral function shows a strong asymmetry about 
the wave vector Ay. In Fig. 14 the spectral function is 
reported as a function of frequency for different wave 
vectors about Ay at T = T C . 



from the Fermi liquid picture appear to be stronger at low 
energy. 

In Fig. 15 the positions of the two peaks of the spectral 
function are reported for different wave vectors about Ay 
at T = T c . The results of our non-self-consistent t-matrix 
approximation (squares and asterisks) are here compared 

with the BCS-like dispersion to = ±^(k) 2 + A^ g (con- 
tinuous and dotted lines), where the BCS gap has been 
replaced by the pseudogap A pg at Ay. It is rather re- 
markable that the coherent peak of the spectral function 
at |k| < Ay gets reflected into the incoherent peak at 
|k| > Ay as the wave vector crosses the "Fermi sur- 
face" (with the characteristic behaviour of an avoided 
level crossing), in such a way that the position of the 
peak at negative frequencies follows almost exactly the 
BCS-like dispersion, provided the value of the pseudo- 
gap is inserted as explained above. A similar behavior is 
displayed by the peak at positive energy. This result sup- 
ports the "heuristic description" of the pseudogap phase 
in terms of an effective BCS approach, that we shall dis- 
cuss in the last Section. 
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FIG. 14. Spectral function at different wave vectors |k| 
about Ay as a function of frequency (in units of £f), for 
{kpap)- 1 = -0.72 and T/T c = 1.001. (Intermediate- to 
weak-coupling regime.) 

It is clear from this figure that for |k| < Ay a well- 
defined peak is found at negative frequencies, and that 
increasing the wave vector to |k| > Ay this peak becomes 
a small and broad incoherent peak. Thus, for |k| < Ay 
the spectral weight of the coherent peak at negative fre- 
quencies decreases as the wave vector k approaches Ay, 
while at the same time the spectral weight of the as- 
sociated incoherent peak located at positive frequencies 
increases, with a transfer of spectral weight from negative 
to positive frequencies upon crossing the "Fermi surface" 
(which is defined as the locus of minimum pseudogap, 
and almost coincides with the sphere k = Ay; note that 
for the coupling value of Fig. 14, Ay is about 10 % smaller 
than A;i?). This clearly shows that the interaction of the 
fermions with pair fluctuations gets increasingly stronger 
upon approaching the "Fermi surface" , so that deviations 
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FIG. 15. Peak positions of the spectral function at nega- 
tive (asterisks) and positive (squares) frequencies versus wave 
vector for (kpap)- 1 = -0.72 and T/T c = 1.001. Full and 
dotted lines represent the BCS-like fit. (Intermediate- to 
weak-coupling regime.) 

To fit the prominent features of .A (Ay , lu) with a sim- 
ple analytic expression (from which the corresponding 
form of X(Ay ,w) replacing Eq. (15) could be extracted), 
we may consider two Lorentians of width 7^ and jr, cen- 
tered at — Al and An, and with weights pl and pn (such 
that pl + Pr = 1 ) , with the labels L and R referring to 
the left and right peaks of A(Ay , u), in the order. In Ta- 
bles I and II we report the values of the fitting parameters 
Ai,Afl,7i, 7ij, and a = 1 — 2pL for the curves of Fig. 13 
(fixed coupling and varying temperature) and of Fig. 16 
(fixed temperature and varying coupling), respectively. 

Note that the asymmetry of the two Lorentians (which 
is controlled by the parameters a and Jl/'Yr) is consid- 
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erable, increasing for increasing temperature or coupling 
(but for the last value of Table II). For temperatures 
and couplings larger than those reported in the tables, 
however, the fit of A(k^ ,uj) with two Lorentians become 
inadequate. 



T/T c 


A L 


Ar 


1L 


7ft 


a 


1.001 
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0.31 
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0.29 
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0.32 


0.28 
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TABLE I. Fitting parameters for the curves of Fig. 13. 
Energy variables are in units of tF- 
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A L 


Ar 


1L 


7ft 
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-1.1 


0.035 


0.035 


0.042 


0.047 


0.060 


-0.72 


0.14 


0.14 


0.11 


0.15 


0.17 


-0.23 


0.47 


0.61 


0.22 


0.33 


0.15 





0.78 


0.84 


0.25 


0.18 


0.008 



TABLE II. Fitting parameters for the curves of Fig. 16. 
Energy variables are in units of 6f- 



Note also that in most cases A^ = = A pg . In 
these cases a relatively simple form for u;) can be 

extracted, yielding: 

£(Ay , oj) = -i(7 - aS) + aA pg 

(l-q 2 )(A 2 g -,5 2 )~2»A pg ,5(l + a 2 ) 

"I t ^~n ~, ft (16) 

oj + aAp g + «(7 + ad) 

where 7 = (7^ + 7l)/2 and S = (7^ — 7z,)/2. Note that, 
even in the symmetric case with a — and <5 = 0, the 
expression ([16]) does not reduce to the form (15) (due to 
the presence of an extra term —ij in Eq. (pq)), unless 
A pg 3> 7 (this condition would be consistent with the 
assumptions under which Eq. (|l^) has been derived only 
when T approaches T c e2i). However, the condition A pg » 
7 is never satisfied by our fits, where A pg and 7 are of 
the same order. 

In Fig. 16 the spectral function at |k| = k^> is reported 
for different values of (kpap)^ 1 from intermediate to 
weak coupling, slightly above the critical temperature. 
Note that, in the weak-coupling regime, the spectral func- 
tion acquires an almost symmetric two-peak structure, 
which differs from the standard BCS result at T = es- 
sentially for the broadening of the peaks due to the finite 
lifetime of the pairs. Note also that the pseudogap near 
the critical temperature decreases with coupling. 

An analysis of the pseudogap opening within a 2D 
attractive Hubbard model in the weak-coupling regime 
has recently been reported in Ref. |3l], by means of the 
non-self-consistent T-matrix approximation formulated 
on the real frequency axis. The frequency dependence 
of the spectral function obtained in that paper (at quar- 
ter filling) resembles the results of our Fig. 16. 

Finally, a comparison of the pseudogap A pg at T c with 
the BCS gap Abcs at T = and with the two-body 
gap eo/2 = €f / (kpap) 2 (which is non-vanishing only for 
Of > 0) is shown in Fig. 17 for all coupling regimes [when 
fi < 0, A BCS (T = 0) is set equal to (^ 2 + A(T = 0) 2 ) 1 / 2 ]. 
Note that in the weak-coupling limit A pg (T — T c ) <C 
Abcs (T = 0) , while in the intermediate-coupling regime 
A pg (T = T c ) ~ A B cs(T = 0). Moreover, in the 
intermediate- to strong-coupling regime (where ap > 0) , 
both Abcs and A pg approach eo/2 from above as the 
coupling is increased. Many-body effects thus increase 
the pair-breaking energy scale with respect to the two- 
body limit. This result resembles the pair-size-shrinking 
effect noticed in Ref. [12] at the mean-field level. 

In this context, it is interesting to mention that, tak- 
ing ep ~ 400 meV as a representative value for cuprate 
superconductors, the range A pg ~ 20 — 120 meV charac- 
teristic of cuprate superconductors corresponds to 0.05 < 
A P g/eF < 0.3, which (as seen from Fig. 17) lies within the 
range identified in Fig. 3 for cuprates. 
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FIG. 16. Spectral lunction at |k| = Ay as a function of 
frequency lj (in units of e F ) for different values of (kpap)^ 1 
and T/T c = 1.001. (Intermediate- to weak-coupling regime.) 
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C. Criterion to distinguish weak from strong 
coupling 

The above systematic study of the single-particle spec- 
tral function from weak to strong coupling suggests the 
following criterion to distinguish by ARPES experiments 
whether a fermion system with an attractive interaction 
lies in the strong- or weak-coupling regime. This criterion 
rests on the analysis of the spectral function at negative 
frequencies (just as determined by ARPES experiments) 
for different values of the wave vector k, and (as dis- 
cussed in the next Section) it is meant to be useful for 
interpreting the experimental data for cuprates in con- 
junction with the two-gap model mentioned in the Intro- 
duction. 

Consider first a system in the intermediate- to weak- 
coupling regime, for temperatures between T* and T c , 



i.e., within the pseudogap region. In this case, the 
chemical potential lies inside the single-particle band 
and almost concides with the Fermi energy. For wave 
vectors smaller than Ay, the spectral function has a 
quasi-particle peak with large spectral weight at nega- 
tive frequencies and a smaller incoherent peak at posi- 
tive frequencies. Upon moving the wave vector across the 
"Fermi surface" (|k| > Ay), the quasi-particle peak shifts 
toward positive frequencies, while the incoherent peak is 
now present at negative frequencies (cf. Fig. 14) and can 
accordingly be measured by ARPES. Restricting to neg- 
ative frequencies and realizing a cut in wave vector space 
which probes the main and the reflected (shadow) bands, 
starting from |k| < Ay ARPES should initially find a 
well-defined quasi-particle peak which, upon increasing 
the wave vector to |k| > Ay, should be reflected as a 
small and broad incoherent peak. Moreover, at |k| = Ay 
the spectral weight at zero frequency remains a sizeable 
fraction of the peak maximum. 

Consider then a system in the intermediate- to strong- 
coupling regime (when jj, < 0), for temperatures between 
Tq and T c . In this case, the chemical potential lies out- 
side the single-particle band. For any wave vector, the 
spectral function has now a quasi-particle peak with large 
spectral weight at positive frequencies and a weaker in- 
coherent peak at negative frequencies. For this reason, 
no appreciable difference in the shape of the spectral 
function should be detected by varying the wave vec- 
tor. Thus, starting, e.g., from k = (k x < 0,0,0) ARPES 
should find a broad incoherent peak which, upon increas- 
ing the wave vectors to (k x > 0,0,0), should not change 
appreciably. In addition, the spectral weight vanishes or 
is much less than the maximum of the incoherent peak 
in a range of frequencies of the order of the pseudogap 
(cf. Fig.9). 

By this token, it is clear that, for a fermionic system 
with an attractive interaction, the wave-vector depen- 
dence and the line-shape of the spectral function at neg- 
ative frequencies have well-pronounced qualitative dif- 
ferences depending on the coupling strength, differences 
which may be detected by a detailed ARPES analysis of 
the spectral function, as discussed next. Recall, however, 
that comparison of our results with ARPES data relies 
essentially on the two-gap model mentioned in the In- 
troduction, and can be complicated by the presence of 
additional sources of quasiparticle scattering in cuprates 
as well as by the fact that the continuum model relates 
strong coupling to low density. Yet, our analysis can be 
useful to understand the evolution of the spectral prop- 
erties along the Fermi surface. 



V. COMPARISON WITH ARPES SPECTRAL 
FUNCTION 

The systematic theoretical analysis of the spectral func- 
tion from weak to strong coupling presented in this pa- 
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per can be used to analyze the spectral intensities mea- 
sured by ARPES in Bi-based superconducting cuprates, 
for which a systematic experimental analysis is also avail- 
able. 

In particular, we consider ARPES intensities measured 
in Bi2212 near the M points of the Brillouin zone as 
well as along the Fermi surface, moving from the M 
points toward the N (nodal) points, in different dop- 
ing regimes and at different temperatures. According 
to our understanding of the phenomenology of cuprate 
superconductors, the effective coupling between fermions 
should increase from the weak- to strong-coupling regime, 
when the doping is reduced from over to optimum dop- 
ing and eventually to underdoping. Moreover, as dis- 
cussed in the Introduction, when moving from N toward 
M points along the Fermi surface, a continuous crossover 
from weakly to strongly coupled fermionic states should 
be observed even at fixed doping. 

We summarize as follows the main results extracted 
from our systematic work in the different coupling 
regimes, which can be compared with ARPES experi- 
ments performed in Bi2212 materials. 

Strong- to intermediate- coupling regime (about M 
points): 

In the strong- to intermediate-coupling regime, where 
the chemical potential is below the bottom of the single- 
particle band (fi < 0), our results show that the spectral 
function is strongly asymmetric in frequency, with two 
peaks, one completely incoherent at negative frequencies 
and the other one coherent at positive frequencies. In this 
case, the wave vectors are meant to be reckoned with 
respect to (one of) the M points. In this regime, the 
prominent features to be compared with experiments are: 

(i) The line shape of the spectral function at negative 
frequencies is quite broad, and the height of the incoher- 
ent peak noticeably decreases with increasing tempera- 
ture (see Fig. 8) or increasing coupling (see Fig. 9). These 
features are in qualitative agreement with the behavior of 
the spectral intensity observed by ARPES experiments in 
the pseudogap phase of underdoped cuprates, by decreas- 
ing doping and increasing temperature. Several ARPES 
measurements show, in fact, that the height of the peak 
in the spectral intensities collected about the M points 
decreases with underdoping, with heavily underdoped 
cuprates displaying a very broad structure with no de- 
tectable peak (cf., e.g., Fig. 2 (left panel) of Ref. |4l] and 
Fig.l (panel a) of Ref. ^, for the doping dependence of 
the spectral weight about the M points). ARPES mea- 
surements for the temperature dependence of the (quite 
broad) spectral intensities about the M points further 
indicate that the spectra are (slightly) suppressed for in- 
creasing temperature (cf. Fig. 2 (panel b) of Ref. ||). 

(ii) The spectral weight near zero frequency is strongly 
suppressed and a real gap opens in the spectral function 
in the strong-coupling regime (see Figs. 8 and 9). Exper- 
imental evidence for a strong suppression of the spectral 



weight near zero frequency can indeed be found, e.g., 
in Fig. 2 (left panel) of Ref. [ll] for (heavily) underdoped 
samples with T c = 56K. 

Intermediate- to weak- coupling regime (between M and 
N points): 

In the intermediate- to weak-coupling regime, the 
chemical potential lies within the single-particle band 
(/i > 0) and the wave vectors are referred to the center 
of the Brillouin zone. In this case, the salient features of 
our calculations to be compared with experiments are: 

(i) A single quasi-particle peak is present in the 
spectral function above the crossover temperature T* 
(see Fig. 13), implying a well-defined Fermi surface. 
Under-, optimally and over- doped cuprates for wave vec- 
tors near the nodal points unquestionably display quasi- 
particle peaks in the ARPES spectral intensities (cf., e.g., 
Fig.l of Ref. H). 

(ii) Approaching the critical temperature from above, 
the interaction between fermions and (damped) pair fluc- 
tuations determines a suppression of spectral weight near 
zero frequency and therefore the opening of a pseudogap, 
characterized by a finite spectral weight at zero frequency 
(see Fig. 13). In addition, the quasi-particle peak dis- 
perses as a function of the wave vector and, as the wave 
vector moves across the Fermi surface, is reflected as an 
incoherent broad peak (see Fig. 14). ARPES spectral in- 
tensities in underdoped cuprates, measured about the N 
points for temperatures between T* and T c , display this 
feature, even though the reflection cannot be accurately 
identified (probably owing to the low spectral weight of 
the incoherent peak). In particular, a spectral weight 
suppression at low frequencies and a finite spectral weight 
at zero frequency has been found by several ARPES mea- 
surements (cf., e.g., Fig.l (panel b) and Fig. 3 (panel a) 
of Ref. |4|). Experimental evidence for the reflection of 
the quasi-particle peak into an incoherent peak has also 
been found by ARPES measurements of the peak along 
the MY direction in the pseudogap phase of slightly un- 
derdoped cuprates (cf. Fig. 2 (panel b) of Ref. 18), for 
which the intermediate- to weak-coupling regime should 
apply. 

(Hi) Increasing the coupling from the weak- to the 
intermediate-coupling regime, the pseudogap evaluated 
at T c increases and the ratio between the pseudogap at 
T c and the BCS gap evaluated at T = also increases 
(see Fig. 17), about coinciding in the intermediate-to- 
strong coupling region. In all underdoped cuprates, and 
for any wave vector, the experimentally determined pseu- 
dogap at T c clearly increases with decreasing doping, and 
in heavily underdoped cuprates it almost coincides with 
the superconducting gap measured at zero temperature 
(cf., e.g., Fig.3 (panel b) of Ref. |). 



16 



VI. DISCUSSION AND CONCLUSIONS 

In this paper, the evolution (from superconducting 
fluctuations to the bosonic limit) of the pseudogap open- 
ing and the spectral function has been studied in a sys- 
tematic way. A system of fermions in a three-dimensional 
continuum, mutually interacting via an attractive con- 
tact potential, has been examined. In this way, the nu- 
merical calculation of the single-particle Green's function 
has been considerably simplified, yet preserving the main 
physical effects underlying the pseudogap opening. The 
pair-fluctuation propagator, the (one-loop) self-energy, 
and the spectral function have been evaluated as func- 
tions of coupling strength and temperature, from weak 
to strong coupling, and analytic and numerical results 
have been presented. 

In the strong- coupling regime, the pair- fluctuation 
propagator has been shown to have bosonic character 
and the line-shape of the incoherent peak of the spec- 
tral function to be strongly asymmetric about its max- 
imum, with its spectral weight decreasing by increasing 
coupling (or decreasing density) and increasing temper- 
ature. In this regime, two crossover temperatures T£ (at 
which the two peaks in the spectral function merge in 
just one peak) and Tq (at which the maximum of the 
incoherent peak crosses zero frequency) have been iden- 
tified, with Tj* > T * > T c and with T * of the order of the 
binding energy of preformed pairs (ARPES experiments, 
however, can only measure Tq). 

In the intermediate- coupling regime, the line shape of 
the spectral function about the "Fermi surface" resem- 
bles the line shape of the spectral intensity (which is, 
in turn, related to the spectral function) measured by 
ARPES in underdoped cuprates between T c and T* for 
different wave vectors. In particular, we have reproduced 
the main features characterizing the ARPES pseudogap, 
namely, a finite spectral intensity at zero frequency and a 
finite pseudogap at T = T c which is of the same order of 
the superconducting gap at zero temperature. We have 
also found that in the intermediate- to weak-coupling 
regime pseudogap effects are present only in a narrow 
temperature range above the critical temperature, a re- 
sult related with the 3D character of the pair fluctuations 
(in 2D this temperature range should, in fact, be consid- 
erably wider). To obtain a wider temperature range for 
pseudogap effects in 3D it is thus necessary to increase 
the coupling as to reach the strong-coupling regime, at 
the price of destroying the Fermi surface. 

In the weak- coupling regime, the pair fluctuation prop- 
agator acquires the diffusive Ginzburg-Landau character 
and the line shape of the spectral function gets progres- 
sively more symmetric as the coupling is decreased. In 
this regime, the two crossover temperatures T{* and T * 
coincide and are of the order of T c , with the pseudo- 
gap closing and filling-in quickly as the temperature is 
increased above T c . 

It is thus clear that the pseudogap already occurs in 



the one-loop approximation for the self-energy, namely, 
the non-self-consistent t-matrix approximation which we 
have adopted in this paper. A similar non-self-consistent 
(as well as a self-consistent) calculation for the spectral 
function has been reported in Ref. [l5|. However, detailed 
comparison of our results with the results of Ref. |l5| ap- 
pears not to be possible, since in Ref. |l5| the shift of the 
chemical potential was not properly taken into account 
when evolving from the weak- to intermediate-coupling 
regimes. 

Most significantly, the results presented in this paper, 
concerning the temperature and wave- vector dependence 
of the spectral function in the pseudogap phase, are in 
qualitative agreement with Monte Carlo simulations of 
the 2D attractive (s-wave) Hubbard model. In particu- 
lar, in Refs. |24| and ^5] the spectral function obtained by 
Monte Carlo simulations is reported in the intermediate- 
coupling regime for different temperatures and wave vec- 
tors. Monte Carlo simulations clearly show that in the 
pseudogap phase the spectral function has a two-peak 
structure, with the incoherent peak smoothly emerging 
from the main peak as the temperature is lowered below 
T* . In addition, moving the wave vector across the Fermi 
surface, the main peak is reflected in a shadow incoher- 
ent peak, as reported in Ref. ^3|. Monte Carlo simulations 
on the 2D attractive (s-wave) Hubbard model thus give 
further support to our non-self-consistent t-matrix ap- 
proximation, suggesting that dimensionality and lattice 
effects do not modify appreciably the main qualitative 
features of the pseudogap phase, obtained by our work 
for a 3D continuum with a contact potential. 

The evolution of the spectral function from weak to 
strong coupling in the pseudogap region can be ratio- 
nalized in a heuristic way by a BCS-like approach at 
zero temperature. The BCS spectral function is given 
by A(k,uj) = u(k) 2 6(uj - E(k)) + v(k) 2 S(io + E(k)), 
where E(k) = ^(k) 2 + A 2 is the quasi-particle energy, 
u(k) 2 = I(l + £(k)/£(k)) and^k) 2 = i(l-£(k)/£(k)) 
are the weights of the quasi-particle poles at positive and 
negative frequencies, respectively. In Fig. 18 we sketch 
this spectral function vs. frequency for three charac- 
teristic wave vectors, namely, |k| >> kp (top panel), 
|k| = kp (central panel), and |k| << kp (bottom panel). 
This behavior of the BCS spectral function for different 
wave vectors is seen to reproduce the qualitative features 
of the spectral function which we have obtained for the 
pseudogap phase above T c in the intermediate- and weak- 
coupling regimes (cf. Figs. 13-16), provided the BCS gap 
A at T = is suitably replaced by the pseudogap A pg 
(obtained from the numerical calculation above T c ) and 
the Fermi energy is replaced by the renormalized chemi- 
cal potential /j,'(T, kpap), which is coupling and temper- 
ature dependent (and also provided the delta-like peaks 
of Fig. 18 are suitably broadened). 
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FIG. 18. Schematic representation of the BCS spectral 
function as a function of frequency for three different wave 
vectors: |k| >> fcp (top panel), |k| = kF (central panel), and 
|k| << kF (bottom panel). 

Note further that even the spectral function in the 
strong-coupling regime can be heuristically understood 
by comparing the BCS spectral function of Fig. 18 for 
|k| 3> kp with our results of Figs. 8 and 9. The dis- 
tribution of the spectral weight and the position of the 
peaks in the spectral function above T c appear thus to be 
in generic correspondence with the features of the BCS 
spectral function at T — 0. The intimate reason why 
the BCS approach at T = is able to give information 
on the pseudogap state above T c remains to be clarified, 
for instance, by extending the present non-self-consistent 
t-matrix from above to below T c . 

Other kinds.|-o£ fluctuation propagators (like, charge- 
density waveE3c3, spin-density waveoca, and phase,fluc- 
tuations above the Kosterlitz-Thouless transitioned) re- 
sult into peak structures in the two-particle Green's func- 
tion and into an associated pseudogap opening in the 
single-particle spectral function. In particular, the pio- 
neering work by Kampf and Schrieffenfa considering anti- 
ferromagnetic fluctuations coupled to fermions has shown 
that the associated spectral function evolves from one 
peak in the Fermi-liquid regime to two peaks in the fluc- 
tuation regime. In addition, the antisymmetric struc- 
ture of the imaginary part of the susceptibility used by 
Kampf and Schrieffer is reminescent of the behavior of 
the imaginary part of our pair-fluctuation propagator in 
the weak-coupling regime only. 

Further detailed ARPES (and, possibly, inverse pho- 
toemission) experiments are awaited to ultimatly distin- 
guish the microscopic origin of the pseudogap in under- 
doped cuprates and to unambiguosly identify the charac- 
teristic features of the spectral function obtained by our 
analysis in different doping and coupling regimes. 
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